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Abstract 

In this work, we calculate the form factors for ,J/^p — > D'*-''^ induced by the flavor changing neutral currents 
(FCNC) in terms of the QCD sum rules. Making use of these form factors, we further calculate the branching 
fractions of semileptonic decays .J/ip D''*^^l^l^ {I ~ e,/i). In particular, we formulate the matrix element 
{J/ilj\Tfj_^\D*'^) with Tfj^iy being a tensor current, which was not fully discussed in previous literature. Our analysis 
indicates that if only the standard model (SM) applies, the production of single charmed mesons at the present 
electron-positron colliders is too small to be observed even the resonance effects arc included, therefore if an 
anomalous production rates are observed, it would be a hint of new physics beyond SM. Even though the predicted 
branching ratios are beyond the reach of present facilities which can be seen from a rough order estimate, the more 
accurate formulation of the three point correlation function derived in this work has theoretical significance and 
the technique can also be applied to other places. In analog to some complicated theoretical derivations which do 
not have immediate phenomenological application yet, if the future experiments can provide sufficient luminosity 
and accuracy, the results would be helpful. 

PACS numbers: 13.20.Gd, 14.65. Dw, 11.55.Hx 



I. INTRODUCTION 



It is widely considered that rare decays of JM can offer an ideal opportunity to study non-perturbative 



QCD effects and the underlying dynamics [11,1^. On other aspect, compared with the extensive studies 
of strong and electromagnetic decays of J/V') both experimental and theoretical investigations of weak 
decays of J /if: are much behind due to their small fractions. 

Thanks to the progress in accelerator and detector techniques, interest in the weak decays of charmo- 
nium is being resurgent. With incomparably large database on J /ip and other ■0— family members, the 
BES collaboration will measure some rare decays including the semi-leptonic [3!] and non-leptonic modes 
J] with high accuracy, and more further theoretical and experimental studies would follow. Theoretically, 
the semileptonic decays of J induced by the flavor changing currents were analyzed in our earlier work 
[^, where the QCD sum rules (QCDSR) approach 



l.0.fl 



a, 



10(] was employed to compute the transition 



form factors. Subsequently, by utilizing the form factors obtained in terms of QCDSR we carried out 



computations on the rates of non-leptonic decays of J /ip [U| under the factorization assumption. Very 



recently, weak decay of J/V' into the final states involving a pseudoscalar meson were also studied by 



authors of Ref. [1^ where the covariant light front quark model was employed, thus their result can be 
regarded as a cross check of that estimated in QCDSR. 

At the quark level, the decay of J /ip induced by the flavor-changing neutral current (FCNC) is realized 
via c ^ u transition, which should be very small due to the Glashow-Iliopoulos-Maiani (GIM) mecha- 
nism [l^, whereas at the hadron level the long-distance effects may have the same order of magnitude. 
Although the FCNC processes for B and K cases are comprehensively studied, the FCNC decays in 



the charmed mesons has not caught enough attention due to the 



str onger G IM suppression for up-type 
[20I . As aforementioned, 



16|, 



13, 



quarks, which is also responsible for smallness of — mixing 
the progress of detection techniques and facilities allows much more accurate measurements on the rare 
decays, so theoretically we need to calculate the production rates and see if the expected precision is 
indeed possible to observe a non-zero fraction at the updated facilities. Thus, in this work, we would like 
to take a step forward to investigate the FCNC processes J/ip ^ D'^*^^1'^1~ in the standard model. 



Following the procedure given in Ref. , we will employ the three-point QCDSR to derive 

the form factors. The QCDSR has been proved to be an effective tool to calculate various hadronic 
matrix elements where non-perturbative QCD effects dominate. The sum rule technique for three-point 
correlation functions was first used to describe the pion electromagnetic form factor at intermediate 



momentum transfer 



22 



23] and hence this approach has been applied to various weak decays 



24 



An alternative approach is the light-cone sum rules where the light-cone distribution amplitudes of 



2 



hadrons are employed 26|, |27| . |28| . |29[ to calculate the form factors in similar processes. In this work, 
we only concentrate ourselves in the QCDSR. To evaluate a transition process, calculation of three- 
point correlation function is needed and obviously it is much more complicated than the calculations of 
two-point correlations. 

The structure of this paper is organized as follows: After this introduction, we will firstly display the 
effective Hamiltonian relevant to the semileptonic decay J /ip ^ and then derive the sum rules for 

the form factors in section [III The Wilson coefficients of various operators contributing to the correlation 
functions are calculated in much detail in section IIIII making use of the operator product expansion 
technique. In particular, the Wilson coefficients of gluon condensate and quark gluon mixing operator 
are dealt with in the fixed-point gauge, i.e., Xfj_A^ = 0. Furthermore, the inputs for the numerical 
computations of form factors are presented at the beginning of section ITVl and then an extensive analysis 
of sum rules of the form factors are performed. We explicitly show the Borel platform where the form 
factors are stable with respect to variations of the Borel masses Mi and M2. The rates of semileptonic 
decays J/ip to D^*^^ are numerically evaluated in section |Vl and the last section is devoted to our 
discussions and conclusions. 



II. THE STANDARD PROCEDURE 

A. Effective Hamiltonian for semileptonic decays of J/t/j to 

The quark level FCNC transition c — > ul~^l~ for the semileptonic decay of J/ip D^*^'^l^l~ is 
described by the effective Hamiltonian 

-2m,C'^ff{^j)uia^,^{l + 75)0/^^], (1) 

with q being the momentum of the lepton pair. In Eq. ([1]), the Wilson coefficients contain the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements. The explicit forms of Cj{/{fi) and Cio can be found in 



literature 



3, 



31 



32 



33 



3j] which are displayed as follows 



C'/f{fl) = Cgifl) + Cn^g, S, /i) + C^g^%Zg, S, /x), 

C'/-f{fi) = Crifi) + Cr{z„ s, 11) + s, ix), (2) 

where the functions Cf°'^{zq, s, ^) and Cl'^^{zq, s, fi) represent the contributions from the continuum and 
resonance parts of self-energy loops of dd, ss and bb and Zg and s are defined as Zg = mg/rric, s = /m^ 



3 



with the subscript q denoting d, s and h quarks. C^°^{zq, s, fj,) caused by the leading order mixing between 



Oi with O9 is given as 



(3) 



q=d,s,b 



with 



2arccot(A/-:2. _ 1) 



In 



ZTT, 



for s < Azl\ 



for s> Azq, 



(4) 



where \q is the CKM matrix element \q = V*qVuq- The contributions of resonances from quark loops 



to Cq^^ ill) can be expressed by Cg'^'^(zq, s, /x) as a shift of the Wilson coefficient C^{ii). Cg'^^{zq,s,iJ,) is 



given m 



37^2 

Cg^^{z, s) = ^ '^i':^ 



mViVvi^i+i- 



(5) 



"v, - q- imv.Vvi 

where m is a free parameter to compensate the deviation caused by the approximation of native factoriza- 



tion 



36 



37 



38l |. and can be adjusted to reproduce the branching ratio of non-leptonic decays D — > ViX 



The numbers of for light vector mesons were calculated in Ref. 
CY°'^{zq, s, fi) and Cj^^{zq, s, ^i) are small and can be neglected. 

For readers' convenience, we collect the Wilson coefficients at = mw as 



32l | as Kp = 0.7, = 3.1 and = 3.6. 



C-j{mw) 
C9{mw) 
Cw{mw) 



q=d,s,b 

^ ^ XqiiC'^'^iXq) + C^iXq)) - 2sm^ Gw (FliXq) + (Xq))], 



sin^e 



w 



q=d,s,b 



1 



sin 9iy 



(6) 



q=d,s,b 



where Xq = mg/m^, Qw is the weak mixing angle. The explicit expressions for Fi{xq), F2{xq),C''°^{xq 



and C^{xq) can be found in Ref. 



4l| and are also included in our Appendix Rl 



In order to obtain the decay rates of J/ip ^ D^*^^, we need to calculate the hadronic matrix elements 



which are usually parameterized in the following forms 
{D'>{p2)\qa^uq''{l+l5)c\J/He,pi)) 



43 



4|: 



2ieppo.pePpfp^Ti{q') - [e^(m^ - mjj) - (e • q)ipi + P2)^]T2{q') 



+ {^■q)[q, 



M 2 2 



{Pi+P2)^mq^), 



(7) 



(5*0(62, P2)|g^T;..g"(l+75)c|J/V(ei, Pi)) 



i{m^ + mz).)(ei • (^*2)[Pif, f 2 oP'^ulT^iq 

(ei • 9) (£2 • q)\Pl^^ — — 2 _ 2 P2^']^5(g 



- i{m^ + mD*)[ei^{e*2 ■ q) - e2^(ei • q)]f(i{q'^) 

+ - - - 2(ei • '7)P2^](e; • q)f7{q^), (8) 

where the totally anti-symmetric tensor is defined as Tr[7^7,y7p7o-75] = ^ie^^pa as a convention adopted 
in our work. It is worth emphasizing that the parametrization of hadronic matrix elements for J/ij: 
decays to vector charmed meson, shown in Eq. dS]) is new and has not ever emerged before. Besides, the 
above parametrization forms are also chosen by the requirement that the stable platform with two Borel 
variables can be achieved to assure our predictions credible. 



B. Sum rules for transition form factors 

1. Sum rules for transition form factors of J /ij: 

As for the FCNC process J/ip ^ D^, both the 'V — ^" current and the tensor operator can contribute 
to the decay amplitude. Here the former one can be directly obtained from the case of J/ip to D^^ by 
exchanging s or d quark into u quark, however, the latter one has not appeared ever before, hence 
we should re-derive the sum rules for the form factors involved in the hadronic matrix element where 
the tensor operator is sandwiched between J/^p and states. Following the standard procedure, the 
three-point function is set as 

where the current ju^^{y) = c{y)^yc{y) represents J/ip channel; j/x(0) = uap,u{l + ^^)q'^c describes the 
weak current for J/ijj to and j^'^ {x) = c{x)i"f^u{x) denotes the Z)^ channel. Inserting two complete 
sets of states with the quantum numbers of J/iJj and mesons simultaneously into the above correlation 
function, one can arrive at the hadronic representation of the three-point function as 

^ ^ fDom3o{D^\ji.\J/^)mj/^fj/^ef ^ r r Ppui^i, S2,q^) 

{m^j/^-Pi){iT^%-vl){'mc + rnu) J Jy.,^ ^ ^ (si - Pi)^(s2 - ^2)^ 
+subtraction terms. (10) 



5 



The subtraction terms are polynomials of either pi or p2, which will disappear after performing the double 
Borel transformation Bj^j2B]yj2, with 

Applying the operator product expansion technique to the in the deep Euclidean region, we achieve 
the expression of this correlation function as 

n^i. = ifoef,ua(3PlP2 + flPl/^Plu + f2P2^P2y + hP2^,Pli, + fiPl^jP2y + h9liu, (12) 

with each coefficient fi contributed from both perturbative part and non-perturbative condensate, i.e., 

/, = /r*I + friqq) + fF^'iGG) + ff'iqGq) + .... (13) 

Comparing the two different expressions for Ilf^y calculated in the QCD and hadronic representations 
and performing the double Borel transformation on variables pi and p2 , we can extract the sum rules for 
the form factors involved in the decay mode of J/^p to as 

Ti(,2) = ^f'- e<l^'h<-"''^MlMlBh, (14) 

Uq') = r-2 T;^?, -e'-y^'^^e'-'r^J^'^MlMlBh, (15) 

Uq') = + e-'V^'^e-'h/^^M^MlBih - h). (16) 



*0 



2. Sum rules for transition form factors of J /ijj D 

Now we are ready to derive the sum rules for the form factors which are responsible for the decay 
channel oi J/ip — > D*^. Now the three-point function can be written as 

n>.up = i' f d''xd^ye-'P^-y+'P---{0\jf'\x)j,{0)j;^/^{ym, (17) 



where the current °(x) = c{x)^pu{x) describes the D*^ channel, and j'u'^{y), J^(0) are the same as 
that in last subsection. The matrix element defined by the "V-A" operator can be gained directly from 
the decay of J /ip to presented in the previous subsection. On the one hand, one can write the 
phenomenological representation of II^i^p at the hadron level as 

m^,ofD*oe'p^' {D*^\jf,\J/'4j)mji^fj/^ef r r p^^^p{si, S2,q^) 



^^'^'^ ("ij/^ - Pi)("^5.o - pi) ^ J Jei2 ~ ^i)(*2 - pI) 

-l-subtraction terms. (18) 



On the other hand, the correlation function Tl^^p can be calculated at the quark level as 

= Fief,uaf3PlP2Plp + F2ef,uapPlP2P2p + Fse^pappfp^Pl^ + heppaf^pfp2P2,y + pcffPlI^Pl p 

+F6e„pai3PlP2P2p + iFrg^uPip + iFsQppPiu + iFggupPip + iFiogpuP2p + iFiigppP2u + iFi2gupP2p 
+iFi'iPipPivPip + iFiiP2pP2vPip + iFi5PipP2uPip + iFieP2p.PiuPip + iFi7P2pP2yP2p 

+iFiSPlpPluP2p + iFi9P2pPluP2p + iF20PlpP2uPlp, (19) 

where each of the above coefficients Fi receives both perturbative and non-perturbative contributions 

Fi = *I + F;"^{qq) + i^«^(GG) + Pf^qGq) + .... (20) 

Finally, equating the above quark-level and hadron- level forms of Hpup, we obtain the sum rules of the 
form factors as 

f^^q^) = - "^l'^ e<'^"e'^l*o""'"MlM^B{Fi - F,), (22) 

(r - "^v- + ^D*o)'^^U'^D*ofD*o 
xB[{Fi + Fs)*?' + {h - h){ml„ - ml)], (23) 

n{q^) = -i— —e<''^'h'^l*o"'''MlM^m, (24) 

f,[q^) = "^£»o - e"'l/''?e"'l*o/''"M?MlBF,,, (25) 

n{q^) = -^-^ — e'^'/^^'e'^i^o/^'^'MfMlMi, (26) 

f,{q^) = , . , ^5*0 - e"^l/^h'-U/''^MfMll3{F, - Fn). (27) 

Now we have achieved the sum rules for the form factors, the next step is to calculate the Wilson coef- 
ficients corresponding to the various operators in the operator product expansion at the deep Euclidean 
region (— ^ 0) in next section. 

III. THE CALCULATIONS OF WILSON COEFFICIENTS 

In this section we calculate the Wilson coefficients. To guarantee sufficient theoretical accuracy, 
the correlation functions are required to be expanded up to dimension-5 operators, namely quark-gluon 
mixing condensate. The dimension-6 operators, such as the four quark condensates, are small and further 
suppressed by O(a^), so can be safely neglected in our calculations. 




id) (e) (/) 




{g) ' ih) (^) 




FIG. 1: Graphs for the Wilson coefficients in the operator product expansion of the correlation function, (a) is for 
the contribution of unit operator; (b) for the two-quark condensate; (c-h) describe the contributions from giuon 
condensate, (i-j) is for the quark-gluon mixing condensate. 



A. Wilson coefficients of the correlation function XI^j^ 

The diagrams which depict the contributions from the perturbative part and nonperturbative con- 
densates are shown in Fig. [TJ The first diagram results in the Wilson coefficient of the unit operator; 
the second diagram is relevant to the contribution of quark condensate, obviously one can neglect the 
heavy-quark condensate at all. The Wilson coefficient of the two-gluon condensate operator is obtained 
from Fig. HJc-h). The last two diagrams in Fig. [Hi-j) stand for the contribution of quark-gluon mixing 
condensate. In this work, all of the Wilson coefficients are calculated up to the lowest order in the running 
coupling constant of strong interaction. 



8 



1. The calculations of perturbative contributions to H^i, 



The contribution of perturbative part to the three-point correlation function 11^,^ comes from Fig. [T] 
(a), which can be expressed as 

Again, we need to express Cj^fJ^ in the form of dispersion integrals. Then, we arrive at the following 
expression 

(si -Pl)(s2 - pi)' 

The integration region is determined by the following condition 

_ ^ < 2s,{s,+ml-ml)-s,{s,+S2-q') ^ ^^^^ 
~ Ai/2(si, S2, g^)A^/2(m2, si, m^) ~ ' 

where X{a,b,c) = a? + b"^ + (? — 2ab — 2ac — 26c. Then, following the standard approach, putting all 
the internal quark lines on their mass shells in terms of the Cutkosky's rules, we can derive the spectral 
density /O^fT* as 

-vert ■ -pert a , -pert , -pert , -pert , -pert , -pert /o-i\ 

Pliu =W0 W/3PlP2+/5l PI11PIU + P2 P2i,P2u + P3 P2i^P1u + Pa PltMP2u + P5 9t^u, (31) 

and the explicit expressions of p^*^^* are collected in Appendix iBl for the concision of the text. 
2. The contribution of gluon condensate to 11^^ 

Now let us focus on the computation of the Wilson coefficient corresponding to gluon condensate. In 
particular, it is worth to emphasize that the contributions of gluon condensate to the correlation function 
no longer vanish, even after performing the double Borel transformation with respect to the variables p\ 
and This point is an important diff'erence between the sum rules of vector current and tensor density. 
The calculations are much the same as that for the case of Il^jy in Ref . , and the only difference is that 
the weak decay vertex "7/^(1 — 75)" is replaced by the tensor one "(T^^g'^(l + 75) ". Besides, we also need 
to rewrite the Wilson coefficient in the form of dispersion integral as that for the perturbative part, i.e.. 



= / / ds,ds2-^^^M^^^^^^, (32) 

^ ' ' {si - pi)[s2 - Pi) 



GG _ / j„ j„_ 



with the integral region being the same as that for the perturbative one. 

The next step is to decompose the above spectral density into various Lorentz structures, namely 

P'ilu = ipo'^etMua0P?P2 + Pl'^Plf.Plu + P2^P2f,P2u + ~Pz^P2p.P\y + ffl^P\p.P2y + '^ff^^, , (33) 



with the exphcit expressions of pp*^ displayed in Appendix [B] for completeness of the paper. The contin- 
uum subtraction should be carried out not only for the perturbative diagram, but also for the contributions 
of the gluon condensate. 



B. Wilson coefficients of the correlation function II^i^p 

Now, we turn our attention to the operator product expansion for the three-point function H^^p in the 
deep Euclidean region, which can be extended to the concerned physical region analytically. Repeating 
the previous procedures but replacing the vertex for a pseudoscalar meson to that for a vector meson, 
one can immediately obtain the expressions of the Wilson coefficients for all the concerned operators. 



1. Calculations of perturbative contributions to Ilp,up 

We can write the perturbative contribution to Ilf_iup shown in Fig. [1] (a) as 

The perturbative part should be expressed in the form of dispersion integral for performing an efficient 
subtraction of the continuum states. In other words, we have 



CZp = ds,ds2 ,'''''\: (35) 



^,.p J J -l-2^,^_^2)(,^_^2)' 

where the integral region is the same as before. Following the standard approach, then, we can analyze 
the spectral function for the perturbative part as below 

-ryprf -.pert \ - pert \ - pert „ g . - pert , - pert 

P]iup = Pi WAP1+P4 WAP2+/'5 ^p-uapPlPlPliy + '^Pr QiiiyPlp + ip'g QvpPlp 

-pert . -^pert 



+ip'n 9ppP2u + 9ypP2^, + (36) 

where only the structures related to the form factors are listed for simplification. Furthermore, the 
explicit forms of /o'f^'^* which are tedious, can be found in Appendix O 



2. The calculation of gluon condensate to 11 



Now we concentrate on the calculations of the Wilson coefficient of gluon condensate for li^up- The 
Wilson coefficient is not equal to zero for the gluon condensate in the operator expansion of Ti^yp. The 
dispersion integral for this Wilson coefficient can be written as 



= / / (37) 

^ ^ < < {si -pi)[s2-pi) 



GG _ j„ f'/^'^P 
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with the integral region being the same as that for the perturbative part. 

Next, we can decompose the above spectral density into various Lorentz structm'es as 



-CO ^ .GG \ ~.GG \ ~.GG ^ /5 , ~ GG . ~.GG 

P/i^p = Pi WAP1+/54 ^fiua/BPlPlPlu + W7 9^lvPlp + ^P9 QvpPl^ 

~ GG ~ GG 

9tipP2u + ip'i5 9upP2p + (38) 
where the explicit forms of p'.^ are given in Appendix [Cl 



IV. NUMERICAL ANALYSIS OF FORM FACTORS IN QCD SUM RULES 



Eventually we are able to calculate the form factors numerically. Firstly, we ex plic itly present all the 



input parameters which are adopted in our numerical computations, as below 



451, 



46 



471 



mdmc) = 1.275 ± 0.015GeV, 
as(lGeV) = 0.517, 
m£,o = 1.865GeV, 



m„(lGeV) = 2.8MeV, 



m 



3.097GeV, 



337ti^MeV, 



foo* = 240tl[jMeV, 



m£,*o = 2.007GeV, 
/do = 166;?oMeV, 
i^C^^C^"') = (0.005 ± 0.004)GeV^ 



(39) 



{uu) ^ -(1.65 ± 0.15) X lO-^GeV^ 



All the QCD parameters are set at the renormalization scale around 1 GeV. To reduce theoretical 
uncertainties in the three-point sum rules of the weak transition form factors, due to masses of quarks. 



threshold parameters and Coulomb-like corrections for J/V' effectively [48l |. we apply the decay constants 



fj/^ and f£)(*)o which are calculated with the two-point QCD sum rules up to the leading order of as, 
to the three-point sum rules. The details about the calculations of the decay constants of both J/ip and 
/){*)o jj^ ii^Q framework of QCD sum rules, are presented in Ref.js]. 

For determining the threshold parameters and demands the QCD sum rules results to 

be sufficiently stable with respect to variation of Mf and M| within relatively large regions, and their 
values should be around the mass squares of the corresponding first excited states. As for the heavy-light 
mesons, the standard value of the threshold in an X channel should be = {fnx + Ax)^, where Ax is 



about 0.6 GeV 



50 



51 



5d, 



53( 1 ■ and we simply take it as (0.6 it 0.1) GeV for the error estimate in our 



50, 



51 



53I ] , we select 



numerical analysis. For the heavy charmonium, following the method given in Ref. 
an effective threshold parameter to ensure the appearance of a satisfactory platform which is around the 
mass square of 'ip{2S). In this way, the contributions from both the excited states including ip{2S) and 
the continuum states are contained in the spectral function. 
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A. Analysis on the sum rules for the form factors 



1. Evaluation of the sum rules for the J/ip — s- 

With all the parameters listed above, we can obtain the numerical values of the form factors. The 
form factors should not depend on the Borel masses Mi and M2 in a complete theory. However, as 
we truncate the operator product expansion up to dimension-5 and keep the perturbative expansion in 
as to the leading order, dependence of the form factors on these two Borel parameters would emerge. 
Therefore, one should look for a region(s) where the results only mildly vary with respect to the Borel 
masses, so that the truncation is reasonable and acceptable. 

With a careful analysis, s? = 13.7 GeV^ and = 6.1 GeV^ are chosen for the sum rules of form 
factors Ti [i = 1,2,3). As commonly understood, the Borel parameters and Af| should not be too 
large in order to ensure that the contributions from the higher excited states and continuum are not too 
significant. On the other hand, the Borel masses also could not be too small for the sake of validity of 
OPE in the deep Euclidean region, since the contributions of higher dimension operators pertain to the 
higher orders in = 1,2). Unlike the treatment adopted in previous literature 



21 



24 ] where the 



ratio of Mi and M2 was fixed, in this paper, when ca 
independently as suggested by the authors of Ref. [4£ 



culating the form factors, we let Mi and M2 vary 



5J. 



As observed in last section, the contributions of gluon condensate are nontrivial for tensor density, 
that is different from the case of vector current. We display the form factors at zero momentum transfer 
in Fig. El As for the form factor Ti, the Borel masses are set as G [6.0, 8.0]GeV^ M| G [1.0, 2.0]GeV^ 
according to the condition that contributions from both the continuum states and the non-perturbative 
gluon condensate to the total sum rules are no more than 40 % and then Ti{q^ = 0) = 0.27l[]:|]^ is resulted 
in. Here we have combined the errors induced by the variations of Borel masses, threshold values, the 
mass of charm quark, decay constants of involved mesons as well as gluon condensate. The value of T2 is 
set as 0.22l|]:g^ by the constraint that contributions of neither the higher states nor the dimension-4 gluon 
condensate can exceed 40 % of the total contribution to the whole sum rules, and it determines the Borel 
region as Ml G [4.0,6.0]GeV2,M| G [1.0, 2.0]GeV2. Finally, the Borel masses for the form factor T3 are 
chosen as G [8.0, 10.0]GeV^,M2 G [1.0, 1.5]GeV^ under the requirement that the contributions from 
both the continuum states and the gluon condensate should be less than 40 %, thus we have T'^{q^ = 0) 
as 0.03T^Q • 

Naturally, we continue to investigate dependence of the form factors at the region q^ G [0, 0.47] GeV^ 
evading the non-Landau- type singularities. The results are shown in Fig. [3] within the given Borel region. 
We fit the form factors with the double-pole approximation for phenomenological applications. Here, 

12 




Mi^(GeV^) 




Mi^(GeV^) 




Mi^(GeV^) 



M2-'(GeV) 



FIG. 2: various form factors Ti,T2 and T3 at = responsible for the decay of J/tp I?"within the Borel 
window. 



one notices that the form factor decreases quickly with the increase of the momentum transfer as 
g2 > O.SGeV^, that is a consequence of the naive and artifact treatment of the continuum density in our 
model and also owing to the smaller gap for the kinematical threshold {q )max = 1.5GeV2 which may 



max 
2 



spoil the operator product expansion [2j]. Besides, we can also find that the q dependence of T3 at the 
region € [0, 0.2]GeV^ is rather mild with the changes of the momentum transfer due to a cancelation 



between the increase of the perturbative part and the decrease of the gluon condensate, the dependence 
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0.1 0.2 0.3 0.4 

q^(GeV^) 



FIG. 3: dependence of form factors Ti, and Tj, for the decay of J/?/" —* D'^ within region without non-Landau- 
type singularities. 



is shown in Fig. [H The other two form factors Ti and T2 can also be written in the double-pole form, 
namely 

^^(^') = 2/ f^L 4/ 4 > (40) 

where the parameters aj and 5j can be determined from the results given by the QCD sum rules in the 
region ^ [0,0.47]GeV2 as 

ar, = 0.75l[J;i, br, = OAOt'oH (41) 



with 



Ti(0) = 0.27tnt T2(0)=0.22tO:[J2. (42) 



2. Estimation of the sum rules for the J/i/j 

Now we numerically evaluate the sum rules for the J/^p D*^ based on the standard method. 
Obviously, the hadronic matrix element for J/ijj D*^ is equal to that for J/i/j D*~ , as long as 
isospin violation effects can be neglected. The threshold used here for the D*^ is also the same as that 
for the D*-, i.e. s§ = 6.8GeV^ To start with, we study form factors at zero momentum transfer with 
the Borel masses presented in Fig. ([SE]). For the form factor Ti, the Borel platform is taken as G 
[6.0,10.0]GeV2,M| G [1.0, 1.6]GeV2 in agreement with the condition that the contributions from both 
the continuum states and the gluon condensate should be less than 35 % of the total contribution, and we 
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FIG. 4: mutual cancelation of dependence of form factor T3 associating with the J/tp ^ decay mode from 
the perturbative part and gluon condensate. 



obtain fi{q^ = 0) as 0.421°;°^. Similarly, we have ^2(^2 = 0) as 0.70l|];[j^ and to be 1.02j:°;j^ with the 
Borel masses being Mf g [6.5, 8.0] GeV^ M| G [1.5, 2.0]GeV2. In addition, the form factor T4 is Q.2Qto^^ 
with the Borel region Mf G [6.0, 10.0]GeV^ M| G [1.5, 2.5]GeV^ Then f^{q^ = 0) with the Borel mass as 
0.41+°:°^ GeV2 and TbCg^ = 0) = 0.38i[!;j]^ with Mf G [6.0, lO.OjGeV^, M| G [1.0, 2.0]GeV2. Ultimately, 
we derive fj{q^ = 0) = 0.1ll[]:[!J with the Borel platform G [6.0, lO.OjGeV^, Af| G [1.5, 2.5]GeV2. 

With the form factors at zero momentum transfer, we can have their values for non-zero . We plot 
the form factors in the kinematical region q^ G [0, 0.42] GeV^ free of non-Landau-type singularities in Fig. 
[71 It can be found that T3 increases quickly as q^ > O.SGeV^, while Ts rises drastically as > 0.2GeV^. 
This point is similar to the behavior of T3 responsible for J/ip ^ , the reason was explained in much 
detail there. The form factors T4 and T^ can be fitted in the single-pole approximation 

(1 - Giq'^/m^Q.o) 

while Ti, and can be written in the following expression 

G^i<l') = n ,2 ' (44) 

(1 - aiq^/m^^,o) 

moreover, Tg is parameterized in the double-pole model 

= 2/ f^°i, 4/ 4 • (45) 

Similarly, the q'^ dependence of T2 is extremely weak, because the dominant contributions of perturbative 
part are almost independent. The parameters and bi can be fixed in terms of the results calculated 
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with QCD sum rules in the region G [0, 0.42]GeV^, then we can extend the above expressions to the 
whole physical region (7^ G [0, 1.2]GeV^. The numbers of these parameters are given as 





- 1 7n+o-i3 

— i.ru_Q2i, 




9 qn+0.06 




T4 






- 2 76+°-28 






- 1 05+0.33 




— Z.i4l_Q j^y, 






— Z.UU_Qog, 






(46) 



and form factors at = are summarized as 



ri(o) 


— U.OO„o.Ol5 


^3(0) 


= 0.801H?, 




^4(0) 


- 16+°-°^ 

— "J-iD-O.Ol' 


^5(0) 


— n Q9+0.02 




^6(0) 






— U.U5y_QQQ3. 


(47) 



B. Discussions on the theoretical uncertainties 

This subsection is devoted to a brief discussion about uncertainties in our calculations. One can 
observe that the errors originating from neither the variations of the Borel masses, nor the thresholds for 
both J /ip and charmed meson channels exceed a level of 10%. The results of form factors are in proportion 
to the inverse of decay constants of the charmed mesons and J /ip, which can be easily observed from the 
definitions of correlation functions. The uncertainties from decay constants of charmed mesons and J/ip 
are at the level of 10 %, which can bring up about 10 % uncertainty to form factors at zero momentum 
transfer, but it would not affect the values of parameters and hi which are used to parameterize the 

dependence of the form factors. Besides, the error bars corresponding to the changes of condensate 
parameters are tiny, since only gluon condensate can contribute to the sum rules, whose effect is at the 
order of a few percents for tensor density transition. Moreover, one can also see that the corrections from 
the light quark masses are not significant because the role of light quark masses is suppressed by a much 
larger energy scale of Borel masses. In addition, the uncertainty caused by the charm quark mass is at 
the one percent level 



45| . We neglect the 0{as) corrections to the perturbative part, which is expected 
to be quite small, and should not result in a drastic shift to the final results. 

In principle, the results of the form factors presented here can be further improved by including the 



non-local quark condensate 21] together with O(^) power correction. The former correction can result 
in a non-vanishing contribution of the diagram (b) in Fig. [1] even after performing the double Borel 
transformation on two variables p\ and p^? however, it almost has no effect on the decay rates, but only 
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FIG. 5: various form factors Ti, T2 and T3 at = responsible for the decay of J/tp D*°within the Borel 
window. 



can moderate the dependence of the form factors. The physical explanation of this effect is that quarks 
in the physical vacuum may have a non- vanishing momentum 2l|] . The latter correction will lead to the 



non-zero contributions of heavy-quark condensate and heavy quark-gluon mixing condensate to the sum 
rules for the form factors. 
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FIG. 6: various form factors T4, T5, Tg and Tj at (7^ = responsible for the decay of J/ip — > £)*''within the Borel 
window. 



18 




FIG. 7: dependence of form factors Ti, T2, T3, T4. T5, Tg and TV for J/V' ^ -D*" within the kinematical region 
without non-Landau- type singularities. 

V. NUMERICAL ANALYSIS OF DECAY RATES FOR J/V' ^ D^*'>°l+l- 

With the form factors obtained above, we can calculate decay rates of semi-leptonic decays J/ip ^ 



D^l + Iv and D'^*^^ll. The related parameters are listed below 



30 
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Gf = 1.166 X 10~5QgV-2, sin^evK = 0.231, 

mw = 80.4GeV, ^^(mfc) = 4.16 ± 0.03GeV, 

me = 0.51 X lO-^GeV, = 0.106GeV, 

^ (48) 

\Vud\ = 0.974, \Vus\ = 0.226 ± 0.001, 

\Vub\ = 3.59 ± 0.16 X 10-3, = 0.226 ± 0.001, 
= 0.973, iKfel = 41.5^1;; x lO'^. 

For the semi-leptonic decay J/V' D^*^^l^l~ {I = e, the differential partial decay rate can be written 
as 

' du, (49) 



dg2 3 (27r)3 32m3 



rmn 



where M = (P5(*)o +Pz+)^ and = {pi+^Pi-)"^] Pi+ andp^- are the momenta of and/" respectively; the 
factor "I" comes from the average of the spin states of J/V'; is the decay amplitude after integrating 
over the angle between the and D^*^^ . 

The upper and lower bounds of u are given as 



19 



Urmn = (i^B(*)o + ^^^f " (V^Bi)o " + ^ E*} - mf^f; (50) 

where )o and are the energies of the charmonium state and the lepton in the rest frame of lepton 
pair respectively 

"2-3, — ^ 



^5(*)o - , - (5ij 

Besides, the explicit form of the decay amplitude M of J/^lJ — > D^*^^l~^l~ is written as 

+ (^(*)0|Cio(/i)^l7M(l - 75)c| J/V')/"7''75/}. (52) 
Then we can obtain the branching ratios for semi-leptonic decay of J/V' — > -D^*] and D^*^^ as 
BR(J/V' ^ ^Oe+e-) = LU^O;!! x lO'^^, BR(J/V ^ ^V/^") = l-OSig:^ x 10-'^ 

(53) 

BR(J/Vj ^ Z)*°e+e-) = 6.30t|^J x 10"l^ BR(J/?/; ^ L>*°^f+^-) = 5.94l|?^ x IQ-^^^ 

where we have combined the various uncertainties for form factors presented in last section into the 
results. As can be observed, the decay rates for the FCNC processes of J/V' ^ L'(*)°/+r should be very 
small, even including the effect from resonances which may enhance the branching ratios considerably, if 
only the SM applies. 



VI. DISCUSSIONS AND CONCLUSIONS 

The weak decays of J /ip meson may serve as a complementary test of the underlying dynamics, 
especially the role of the FCNC in weak decays compared with the strong and electromagnetic processes 
which absolutely dominate the J lifetime, even though it is very difficult to be experimentally observed. 
Due to the progress of detection facilities and techniques, it might be feasible to measure so small 
branching ratios with relatively clean background and huge database in the future. Of course, if the 
measurement is feasible, it would be a better platform for determining the CKM entries because of 
absence of contamination from the spectator. 

Moreover, the rare weak decays of J may be sensitive to the new physics beyond the standard 
model. Once such weak decays were observed with a sizable branching rate in the future colliders, it 
would be a clear signal of new physics effects. 

In this work, we calculate the weak decay rate of J/il) D^*^^l^l~ which is realized via FCNC- 
induced processes in the framework of SM, we find that the rate is too small to be observed in the 
facilities available at present. Namely, our numerical results show that the branching ratios of such 
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decays are at the order of 10~ . In the calculations, we have used the QCD sum rules and taken into 
account possible uncertainties coming from both theoretical and experimental sides. Even though the 
predicted branching ratios are beyond the reach of present facilities which can be seen from a rough 
order estimate, a more accurate formulation of the three point correlation function derived in this work 
has theoretical significance and the technique can also be applied to other places. In analog to some 
complicated theoretical derivations which do not have immediate phenomenological application yet, if 
the future experiments can provide sufficient luminosity and accuracy, the results would be helpful. 
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APPENDIX A: THE DETAILED EXPRESSIONS OF BASIC FUNCTIONS RELATED TO 
THE FLAVOR-CHANGING NEUTRAL CURRENT PROCESSES 



This appendix is devoted to the collection of the basic functions associating with flavor-changing 
neutral current processes, which are taken from 3j] as 



C^ixq) 



1 ^ XqlnXq 



1 K-l)2 



Xg 3 1 



3 2x„ Xq 



4 8xg-l 64 (xg - 1)2 



9' 



Qq 
+ 



1 1 



+ 



13 



Uxg-1 12(a;g-l)2 2{xg - 1)^ 



2 1 



+ 



1 



1 



+ 



1 1 



3Xg-l \3{Xg-iy 6(Xg-l)3 2{Xq-l) 

7 1 13 1 11 

3Xq-l ^ 12 {Xq - 1)2 ~ 2 {Xq - if 

1 1 35 1 5 1 11 



4 1 



Inxr. 



+ 



F2{Xq) = -Q 



6Xq-l 12(Xg-l)2 Q(xq-lf 2 {Xq - l)^ 

1 3 1 



1 1 3 

+ 



AXq-l 4{Xq-lf 2{Xg-iy^ 



3 XqinXq 



+ 



1 1 



+ 



9 



1 



+ 



1 



2Xq-l A{Xq- 1)2 2 {Xq - I 



3 x^lnx 



2 {Xq - 1) 
q — q 



4 {xq - ly 



(Al) 



where Qq is the charge of down quarks q = d,s, b. 
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APPENDIX B: THE WILSON COEFFICIENTS FOR ft^ 



This appendix is devoted to the collection of visible Borel transformed forms of Wilson coefficients 
responsible for the tensor current transition corresponding to the decay of J/V' to as been presented 
in Eq. (jl4N16|) . It can be observed from the text that both perturbative diagram and gluon condensate 
diagrams contribute to the correlation functions non-trivially, which is distinct from the chiral current 
transition remarkably. In the mathematical language, it can be written as 



/, = /r*I + /f ^'(GG) + Oias) + 0(l/m,), 



(Bl) 



where ff^^^, /['^ can connect with p'^^'^ and in light of the following formulae 



-pert 



GG 



fi 



'pert 



or 



fGG 



{mc+rriu)^ 
(mc+mu) 



ds2 I dsi * 



(si -pi){s2-pi) 

ds2 / dsi- 2T7 2T' 



{mc+rriu)' 



J(mc+mu)2 Js^ ^y^i JyJ-2 



The lowest bound of si, i.e., can be determined by the Eq. (j3U|) as 
1 



2mg 



4 /n 2 I , 2\ 2 , 2 I 2 2/ , 2\ 

- (2m + S2 + q + m„ + S2q - m [S2 + q ) 



according to the Landau equation 
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(B2) 
(B3) 

(B4) 
(B5) 



+ \Jm'^ - 2(m2 + S2)m1 + {rrfi^ - 82)'^ ^Jm^ - 2{m'^ + q'^)ml + (m^ - q2)2 ^ 

58 1 . The detailed expressions of pf^'^^^ and pf^ are given as 



(B6) 



-pert I ■z\ 

Po [si,S2,q ) 



-pert I 2\ 

Pi [si,S2,q ) 



:{2sim^ - (4sim^ + {si - 82)^ - A)m^ + 2mumc\ + 2m^si + m^((si - 82)^ - A) 



8^2A3/2 

+ S2)(A - (si - S2f) + g^(2((si + S2)ml + sl + sl + S1S2 - ml{si + S2)) 
-{ml -ml + si+ S2)q^)], 
pP^'\si,S2,q^) = ^^^^{3si(A2 - {2ml " ^^l - si + ^2 + q^fX - 2si{-2ml + 2ml + si - ^2 - q^)X 
+2(si - S2 - g^)(6si?n^ + 2(A - 3si{2ml + si - S2))ml 

.2 , „_^_2 , s2)^)+sig^(6m^-6m^-2si+4s2 + g^ 

l2 o„ r „ I „ „2n I ^^2 ^,2\/„ I „ ^2 



+si(6m„ + 6(si - S2)m-u + (si - S2) ) + sig (6771^ - 6m„ - 2si + 4s2 + )))}, 
^--^^^{3((2mc(m„ - mc) + si)A^ - 2si(s2(-si + S2 - g^) + {ml - m^)(si + S2 - Q'^))A 
-si(si - S2 - q^){-2ml + 2m^ + si - §2 - 9^)A - 2(si - S2 - q^){{si + S2)\ml 
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+si(3(si + S2)m^ - 2(3(si + S2)m^ + (si - S2)(si + 2s2))m^ 
+{si - 32^32 + 3m^(si + 32) + 2m^(si - S2)(si + 2s2)) 

+q'^{-Xml + si(s| + {-2ml + + si)s2 + {ml - ml){-3ml + 3ml + 4si)) 

-2si(m^ -m^ + S2)g^)))}, 

yof''''*(si,S2,g^) = ■^;^Jj^{3{{mc-mu)X{mu3i+mc32-mcq^) 

+ - S2 - q^){\ml + (m^ - ml)3i{ml - - si + S2) + si("i^ - + 82)0'^))}, 
647r^(gi -g^) 
A37^ ' 



p^^(si,S2,g') 



^ac^ 2n 128^2g^((g^_g2)(s^ + 2s2)+g2(-2si + S2 + g^)) 

P2 1^1) j = , 

-GG. 2n 647r2(si(si - S2)(si + 5s2) - g^(3s? + 63231 + si + g2(-3si - 2s2 + 

P4 lSl)'S2,g j = , 

.GG. , 2^ _ l6TT^{3l+32-q^){-3l+32+q'') 

P5 {Sl,32,q ) • (o7j 

APPENDIX C: THE WILSON COEFFICIENTS FOR H^^p 

This appendix is concentrated on Wilson coefficients relevant for the tensor density transition of J/ip 
to decay appeared in Eq. (j21H27p after doing the double Borel transformation. As mentioned before, 
both the perturbative and gluon condensate parts are nonzero in the operator product expansion of 
three-point function accounting for the tensor operator's matrix element, which can be written as 

Fi = *I + F,^^{GG) + 0{as) + ©(l/m^. (CI) 

The connections of F?"^''* and F[^^ with p'f'^''*, p'f^ can be expressed as 



* = / " d32 I 'dsiT^^^-^s^^F^^ (C2) 



-^pert 2 



1 , Pi {si,32,q- 



^;°°=r" ^iT^Sfiif^. (C3) 

J(mc+m„)2 JsL [31 - p(){32 - Pi) 



or 



(C4) 



i3i^^^= r° d32 f'^ d3'Ge-^'/^?'e-^^/'^^lp'f''{ (C5) 

with the lower limit of integrals sf" defined as before. Besides, the obvious forms of p'^^^^ , p'f^ can be 
displayed as 

p'i'"'*(si, 32, q^) = - ^^^^g^^ {(-2(mc + m„)si(si - S2)(m^ - - si + S2)(mc - ruu)^ - rUcX^ 
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-(2(si - 2s2)ml + 2m„(s2 - 2si)m^ - (s? + - 2(m^ + Si)s2)mc + 2m^si)A 
+q^{2{mc - m„)si(m^ - 2(m^ + si - S2)m1 + + 2m^(si - S2) + S2(s2 - si)) 
-2mc{{mc - mu)mu + si + S2)A + (2(mc - m„)si(mc - "i^ + S2) + 'tT'cA)^^))} 
/9'4'"'*(si,S2,g^) = -— ^^^{(4si(si - S2)m^ + 4m^,si(s2 - si)m^ 



+2si(3A - (4m^ + 3si - 3s2)(si - S2))"^|? + 2m„si((4m^ + 3si - 3s2)(si - S2) - A)?n; 
+(A2 - (^2 - 4siS2 + 3si(2m^ + si))A + 2si(si - S2){ml + si - S2){2ml + si - S2))m, 
-2m„si(m^ + si - S2)((si - S2)(2m^ + si - S2) - A) 

+q'^(4(mc - mu)si{m^ + (-2m^ + si + S2)m^ + + Si(s2 - si) - m^(si + S2)) 
+2mc(si + S2)A + (2(mc - m„)si(m^ - + si + S2) - mcA)g^))} 



/f's {si,S2,q) = - S2)((-2mg + 2m„m^ + S2Wtc + muSi)X 



47r2A' 

-2(mc - m„)^(mc + m„)si(m^ - - si + S2)) 

+g^(2(mc - mu)si{m^ - 2{m\ + si - S2)m1 + + 2m^(si - S2) + ^2(52 - si)) 
+ (2m^(mc - m„) - (ttt-c + m„)si)A + (2(mc - m„)si(m^ - + ■52) + "icA)g^))} 
P'7'\si, S2, q^) = - g^^^gy^ {((mc + m„)((mc - m„)^ - S2)A - (rwc + m„)(s2(-si + S2 - 5^) 
+(m^ - + S2- q^))isi - S2 + q^) 

+2K - -.)(Am^ + - ml)s,{ml - ml - s, + S2) + s,{ml -ml + S2)q'))] 
p'l''''\si,S2, q^) = g^^^3y^ {((me - m„)(s2(-si + S2 - 9^) + [ml - mi){si + S2- q^)){si - S2+ q^) 

+\{ml - m^ml - mlm, " + m^{ml + s, - ^2)))}, 
P'u {si,S2,q^) = ^^^jy^^{(mc + m„)si(A + (2(mc - mti)^ + si - S2 - g^)(2m^ - 2m^ - si + S2 + g^))}, 

P'ir*(si> ■52, q^) = - ^^2x5/2 ^i^c>? - {mc - mu)is2{-si +S2- q^) + {ml - m^)(si + S2 - q^))\ 
-{rric + ms){s2i-si + S2 - q^) + {ml - m^)(si + S2 - q^))\ 
+ {mc - mu)si{2ml - 2m^ - si + S2 + g^)A 

-2(mc - mu){{si + S2)Am^ + si(3(si + S2)m^ - 2(3(si + S2)m1 + {si - S2)(si + 2s2))"ic 
+ (si - S2)^S2 + 3m^(si + S2) + 2ml{si - S2)(si + 2s2)) 

+q'^{-Xml + si(s^ + {-2ml + + si)'52 + {ml - ml){-3ml + 3ml + 4si)) 
-2si(m2-?n2 + S2)<?^)))} 
/o'l (si,S2,g^) = -^^{128(mc - m„)7r^(g^ - (2si + 52)17^ + si(si - S2))} 

/o'4 (si,S2,g^) = --^^{256(mc - mti)7r^si(si - S2 + 2g^)} 
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■GG 



P'5 isi,S2,q ) 



A5/2 



~,GG 



P'7 (si,S2,g ) 



A5/2 



{128(mc - m„)7r2(si - S2 - g^)(2si - 53 + g^)} 
{64(me - m„)7r2(si(si - s^f - ^^(Ss? + 12siS2 + + g^(-3si - 2s2 + g^)))} 



•GG 



P'g (si,S2,g ) 



1 



•,GG 



P'll (si,S2,g ) 



A5/2 
1 



•,GG 



P'l5 {si,S2,q ) 



{64(mc - mu)7r^(si - S2 - - S2)(2si + S2) - q^{si - 2s2 + q^))] 

{128(me - m„)^2^i((si - ss)^ + 4g4 - 5(si + S2)q^)} 
{128(mc - m«)7r2((s? + 10siS2 + sl){si - S2f 
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